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Abstract 

We consider a single server retrial queue with the server subject to interrup- 
tions and classical retrial policy for the access from the orbit to the server. We 
^ \ analyze the equilibrium distribution of the system and obtain the generating 

I ' functions of the limiting distribution. 

(N . 

(N ■ 1 Introduction 

o 

o\ . 

Q , Queueing systems with retrials of the attempts are characterized by the fact that an 

arrival customer who finds the server occupied is obliged to join a group of blocked 
customers, called orbit, and reapply after random intervals of time to obtain the 
service. These systems are useful in the stochastic modeling of much of situations 
in practice. We can find them in aviation, where a plane which finds the runway 
occupied remakes its attempt of landing later and we say in this case that it is 
in orbit. In Telephone systems where a telephone subscriber who obtains a busy 
signal repeats the call until the required connection is made. In data processing, 
we find them in protocol of access CSMA/CD. They appear in the modeling of the 
systems of maintenance and the problems of repair among others. For details on 
these models, see the book of Falin and Templeton [4J or the recent book by Artalejo 
and Gomez- Corral [2|. 

We study in this article single server retrial queues with various types of inter- 
ruptions of the server. From a practical point of view, it is more realistic to consider 
queues with repetitions of calls and the server exposed to random interruptions. 
Queueing models with interruptions of service proved to be a useful abstraction in 
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the situations where a server is shared by multiple queues, or when the server is sub- 
ject to breakdowns. Such systems were studied in the literature by many authors. 
Fiems et al. [6] considered an M/G/1 queue with various types of interruptions of 
the server and our work is a generalization to the case of retrial queues. White 
and Christie [T7] were the first to study queues with interruptions of service by con- 
sidering a queueing system with exponentially distributed interruptions. Times of 
interruptions and services generally distributed are considered by Avi-Itzhak and 
Naor ^ and Thiruvengadam [Hj. Other generalisations were considered in the lit- 
erature ( phase-type [5], approximate analysis [15], Markov-modulated environment 
[13] and [in], processor sharing [H]). Gaver [S] considers the case where the service 
is repeated or repeated and begin again after the interruption. 

We consider in this paper a single server retrial queue with server interruptions 
and the classical retrial policy where each customer in orbit conduct its own attempts 
to get served independently of other customers present in the orbit. We can then 
assume that the probability of a retrial during the time interval {t,t + dt), given 
that j customers were in orbit at time t, is jddt + o[dt). Kulkarni and Choi [H] 
studied a single server linear retrial queue with server subject to breakdowns and 
repairs and they obtained the generating functions of the limiting distribution and 
performance characteristics. Artalejo [1] obtained sufficient conditions for ergodicity 
of multiserver retrial queues with breakdowns and a recursive algorithm to compute 
the steady-state probabilities for the M/G/1 linear retrial queue with breakdowns. 
The detailed analysis for reliability of retrial queues with linear retrial policy was 
given by Wang, Cao and Li [16j . 

The remainder of paper is organized as follows. In the following section, we 
describe the model and give the necessary and sufficient conditions so that the system 
is stable. In section 3, we analyze the equilibrium distribution of the system in study. 

2 Model Description 

Consider a single server queueing system in which customers arrive in accordance 
with a Poisson process with arrival rate A. If at the instant of arrival the customer 
finds the server free, it takes its service and leaves the system. Otherwise, if the server 
is busy or in interruption, the arriving customer joins an unlimited queue called orbit 
and makes retrials for getting served after random time intervals. We consider the 
classical policy where each customer in orbit conducts his own attempts to obtain 
service independently from the other customers present in the orbit. We can then 
assume that the probability of a retrial during the time interval {t,t + dt), given that 
j customers were in orbit at time t, is j6dt + o [dt). Service times constitute a series 
of independent and identically distributed (i.i.d.) random variables with common 
distribution function B{t), density function b{t), and corresponding Laplace-Stieltjes 
transform (LST) /3{s) and finite first two moments = (— l)'^ /5*-'^''(0), k = 1,2. 
Interruptions of the service may occur according to a Poisson process with rate u if 
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the server is busy and this type of interruption can be disruptive with probabihty pd 
(or rate = Pd^^) or non- disruptive with probabihty Pn = ^ — Pd (or rate z/„ = p„z/) . 
In the case of a disruptive interruption the customer being served repeats his service 
at the end of the interruption, in the other type the customer continues his stopped 
service. If the server is idle, another type of interruptions may occur according to a 
Poisson process with rate z/j. We call this type idle interruption. The lengths of the 
consecutive disruptive (non- disruptive, idle time) interruptions constitute a series 
of i.i.d. positive random variables with distribution function Bd{t) {Bn{t), Bi{t)), 
density function bd(t) ipnit), bi{t)), corresponding Laplace-Stieltjes transform (LST) 
/3^(s) (/?„(s), /?j(s)) and finite first two moments (Z^^, /3l) k = 1,2. 

Denote by N{t) the number of customers in orbit at time t. Let C{t) be the 
state of the server at time t : C{t) = F if the server is free (and functions normally), 
C{t) = S if the server is busy (and functions normally), C{t) = D if the server is on 
a disruptive interruption, C{t) = N if the server is on a non-disruptive interruption, 
C{t) = I if the server is taking an idle interruption. We introduce the random 
variables ^{t), (.uit), ^^(t) and C/(t) defined as foUows. If C(t) = S then C,(t) 
represents the elapsed service time at time t; if C{t) = D, then ^£)(t) represents 
the elapsed disruptive interruption time at t; if C(t) = N, then ^jv(t) represents 
the elapsed non-disruptive interruption time at t; and if C(t) = I then ij{t) is the 
elapsed idle interruption time at t. 



3 Stability Analysis 

We first study the condition for the system to be stable. The following theorem 
provides the necessary and sufficient stability condition. 

Theorem 1 The system with classical retrial policy and interruptions is stable if 
and only if the following condition is fulfilled 

^ii^^(l + ../3f + .„/3?)<l. (1) 

T^dPiVd) 

Proof. Let {s„; n G N} be the sequence of epochs of service completion time. We 
consider the process F„ = (A^(s„+), C(s„+)) embedded immediately after time s„. 
It is readily to see that G N} is an irreducible aperiodic Markov chain. To 

determine the stability of the system it remains to prove that n G N} is ergodic 
under the suitable stability condition. Let us first consider the generalized service 
time 5* of a customer which includes, in addition to the original service time S of 
the customer, possible interruption times during the service period of the customer. 
Fiems et al. [6] showed that the generalized service time has the Laplace transform 

[S^U- VnPn{s)\ (3 {s + U - I/„/5„(s)) 



[s + u- i/„/3„(s)] - Udf^dis) [1- (3{s + u- z/„/3„(s))] ' 
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hence its expected value is given by 

ES = -PiO) = ^-^—P^ (1 + u,f3i + ^./??) . 

For the sufficiency, we shall use Foster's criterion, which states that a Markov chain 
{Yn] n G N} is ergodic if there exists a nonnegative function f{k), k eN, and 6 > 
such that for all A; 7^ the mean drift 

Xk = E[fiY,,^,)-f{Y,,)\Y,, = k], (2) 

satisfies Xk ^ ~^ E [f{Yn^i) \ Yn = 0] < 00. If we choose f{k) = k we obtain 

E I r„ = 0] = XES = ^(^-J^^d)) (1 ^ ^^f^d ^ ^^^n) ^ 

and we can easily check that 

Xj, = XES-1 = [A (1 - /5(Z/,)) /UdPiUd)] (1 + UdPt + Vn(3l) - 1. 

If we set 

then the condition ([1]) is sufficient for ergodicity. 

To prove that the condition ([1]) is necessary, we use theorem 1 of Sennot et al. [T2] 
which states that if the Markov chain {Y^] n G N} satisfies Kaplan's condition, 
namely Xa; < 00 for all /c > and there is an k^ such that X/c — k > k^, then 
[Yn] n G N} is not ergodic. Indeed, if 

A(1-/?K)) , , „ . 



i^dPii^d) 



(l + z/,/?f + i/„/5^) > 1 



then for f{k) = k, there is a /cq such that pij = for j < i — ko and i > 0, where 
P = (Pij) is the one-step transition matrix associated to n G N}. 
The stability of the system follows from Burke's theorem (see Cooper [7J pl87) since 
the input flow is a Poisson process. ■ 



4 Steady-state analysis 

We investigate in this section the steady-state distribution of the system. Define 
the functions i^o^x), t^N^y) Hii^x) as the conditional completion rates for 

service, disruptive interruption, non-disruptive interruption and idle interruption, 
respectively, i.e., ix{x) = b{x) / {1 - B (x)) , fic>{x) = hd{x) / {I - Bd{x)) , fi^ix) = 
hn{x)/ (1 - Bn{x)) and = hi{x)/ (1 - Bi{x)) . 
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We now introduce the following set of probabilities for j > 0: 





P{N{t) 




C(t) 




PB,j{t,x)dx = 


P{N{t) 


= J, 


C(t) 


^ S, X < ^{t) < x + dx} , 


PD,j{t,x)dx = 


P{N{t) 


= J, 


C(t) 


^ D, X < ^ij{t) <x + dx} , 


PN,j{t,x,y)dy = 


P{Nit) 


= 3, 


Cit) 


^N, m=x, y <CNit) <y + dy}, 


pij{t,x)dx = 


P{N{t) 


= J, 


C{t) 


^ I, X < ^j{t) <x + dx) . 



where i > and y > 0. 

The usual arguments lead to the differential difference equations by letting t 

+00 



+ 00 +00 



d_ 
dx 



{X + j9 + Ui)pFj = J n{x)pB,j{x)dx + j iJ,j{x)pjj{x)dx, (3) 



+ 0O 

+ X + iy + ii{x)^PB,j{x) = {1 - 5oj) XpB,j-iix) + J iiMPN,j{^^y)dy,{4:) 



(^-^ + X + I^d{x)^ PDjix) = {l-6oj)\pD,j-li^), 

(^-^ + X + HN{y)^PN,j{x,y) = {l-SQj)XpN,j-i{x,y), 
■^ + X + fij{x)^Pij{x) = {l-5oj)Xpij-i{x). 



(5) 
(6) 
(7) 



With boundary conditions 



+ 00 

Psj(O) = {j + l)9pF,j+i + XpF,j+ J i^D{x)pD,jix)dx, 



(8) 



Pd,j(0) ^ j PB,j{x)dx, (9) 



PN,j(x,0) = UnPB,j(x), (10) 
Pl,M = ^iPF,j- (11) 

The normalising equation is 

^PF,j + Yl / PB,j(^)dx + Y^ / PD,j{^)dx + Y^ / / PN,j{x,y)dxdy 

3=0 3=0 i 3=0 i j=o i i 

+ Pi,j(^)dx^l 



3=0 
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Define the generating functions 

oo 

j=0 

oc 

Pb(x,z) = ^PB,j{x) z\ 

3=0 

oo 



Pd{x,z) = ^pD,j(x)z^, 

3=0 

oo 

PN{x,y,z) = ^pN,j{x,y) z^, 



3=0 

oo 

Pi{x,z) = ^pij{x)z^, 

3=0 

for \z\ < 1. 

We introduce h{z) = \y — i^nPn ~ '^^) + ^ ~ ^-2] and 
X {z) = h{z) — I'dPd (-^ ~ '^^) (1 ~ (^(-2))) to simplify notation. 
We have the following theorem 

Theorem 2 In steady state, the joint distribution of the server state and queue 
length is given by 

Pf{z) = Pf(1) exp I j ■^{u)dv^ , 

Pb{x,z) = PBiO,z){l- B{x))cxp{-h{z)x} , 

Pd{x,z) = Pb{0, z) ~ (1 - Ba{x)) exp {- (A - A^) x] , 

n[z) 

PN(x,y,z) = PB(0,z)un(l- Bn(y))(l- B{x))exp{-h(z)x}exp[-{X- Xz)y], 
Pi{x,z) = UiPF{z){l-Bi{x))&q>[-{X-Xz)x], 



where 



(1 + iy,/3\) [uli] (ud) + XvnPl (1 - P{yd)f - XiynP>d (1 - Pii^d))] ' 
Xh{z)P {h{z)) - [A + z/, (1 - (3^ (A - A^))] X {z) 



9{zx{z)-h{z)l3{h{z))) 
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